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Abstract 

Th aim of this paper is to study the nonparametric regression esti- 
mators on the sphere built by the needlet block thresholding. The block 
thresholding procedure proposed here follows the method introduced by 
Hall, Kerkyacharian and Picard in [53] , [53] , modified to exploit the prop- 
erties of the spherical standard needlets. Therefore, we will investigate on 
their convergence rates, attaining their adaptive properties over the Besov 
balls. This work is strongly motivated by issues arising in Cosmology and 
Astrophysics, concerning in particular the analysis of cosmic rays. 
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1 Introduction 



Over the last years, wavelet techniques have been used to achieve remarkable 
results in the held of statistics, in particular in the framework of minimax estima- 
tion in nonparametric settings. The pioneering work in this area was provided 
by Donoho et al. in where authors proved that nonlinear wavelet estima- 
tors based on thresholding techniques attain nearly optimal minimax rates, up to 
logarithmic terms, for a large class of unknown density and regression functions. 
Since then, this research area has been deeply investigated and extended - we 
suggest for instance [23] as a textbook reference: specifically, we are focussing on 
the wavelet block thresholding procedure, among the other techniques. Loosely 
speaking, this method keeps or annihilates blocks of wavelet coefficients on each 
given level (for more details, see 23 ), hence representing an intermediate way 
between local and global thresholding, which fix a threshold respectively for 
each coefficent and for all the set of them. This procedure, initially suggested 
in [15] for orthogonal series based estimators and later applied by [24] for both 
wavelet and kernel density estimation on E (see also [IS]), was used in [5J jointly 
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to Oracle inequalities; overlapping block thresholding estimators were studied 
in [6]. The block thresholding was also applied to study adaptivity in density 
estimation in [DJ , a data-driven block thresholding procedure for wavelet regres- 
sion is instead investigated in [7J, while wavelet-based block thresholding rules 
on maxisets are proposed by pQ. 

Even if a huge number of results concerns estimation with the thresholding 
paradigm in standard Euclidean frameworks, such as R or M n , more recent 
applications are being established in more general settings, such as spherical 
data or more general manifolds. In particular, we aim to a highly successful 
construction of a second-generation wavelet system on the sphere, the so-called 
needlets. The needlets were introduced by Narcowich, Petrushev and Ward in 
[36] . [37] ; their stochastic properties, when exploited on spherical random fields, 
were studied in [2] , [3] , [32] and [33] . This approach has been extended to more 
general manifolds by [2DJ , [2T] , [22] , while their generalization to spin fiber bun- 
dles on the sphere were described in [Tg], [ID]. Most of these researches can 
be motivated in view of their applications to Cosmology and Astrophysics: for 
instance, a huge amount of spherical data, concerning the Cosmic Microwave 
Background radiation, are being provided by satellite missions WMAP and 
Planck, see [38], [35], [35], [ID], 00], gT], [ID], 02], [13] and Q3] for more de- 
tails. The applications mentioned here, however, do not concern thresholding 
estimation, but rather they can be related to the study of random fields on the 
sphere, such as angular power spectrum estimation, higher-order spectra, test- 
ing for Gaussianity and isotropy, and several others (see also [8]). As another 
example, we mention experiments concerning incoming directions of Ultra High 
Energy Cosmic Rays, such as the AUGER Observatory (http://www.auger.org). 
The Ultra-High Energy Cosmic Rays are particles with energy above 10 18 eV 
reaching the Earth. Even if they were discovered almost a century ago, their 
origin, their mechanisms of acceleration and propagation are still unknown. As 
described in g], see also [TjJ, an efficient nonparametric estimation of the den- 
sity function of these data would explain the origin of the High Energy Cosmic 
Rays, i.e. if it is uniform, they are generated by cosmological effects, such as 
the decay of the massive particles generated during the Big Bang, or, on the 
other hand, if it is highly non-uniform and, moreover, strongly correlated with 
the local distribution of nearby Galaxies, it means that the Cosmic Rays are 
generated by astrophysical phenomena, as for instance acceleration into Active 
Galactic Nuclei. Massive amount of data in this area are expected to be avail- 
able in the next few years. Also in view of this application, the needlet approach 
was recently applied within the thresholding paradigm to the estimation of the 
directional data: the seminal contribution in this field is due to g] , see also [28] , 
[27] . while applications to astrophysical data is still under way, see for instance 
[IB] . [T7] and [2DJ. Minimax estimators for spherical data, outside the needlets 
approach, were also studied by Kim and coauthors (see [3D], [2DJ, [31]). Further- 
more, adaptive nonparametric regression estimators of spin-functions, based on 
spin pure and mixed needlets defined in [18], [19], were investigated in [12]. In 
this case, the needlet nonparametric regression estimators were built on spin 
fiber bundles on the sphere, i. e. the function to be estimated does not take as 
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its values scalars but algebraic curves living on the tangent plane for each point 
of the sphere. 

This work, hence, extends the results established in [3] and [T2] towards the 
needlet block thresholding procedure following two main directions. First of 
all, we will suggest a construction of blocks of needlet coefficients, exploiting 
the Voronoi cells based on the geodesic distance on the sphere. Then, we will 
define the needlet block thresholding estimator, whose we will achieve an optimal 
convergence rate. In view of this aim, we will use both the needlet properties 
established in [3B], [37J (see also [33]) an d a set of well consolidated standard 
techniques, introduced by [IT] (see also [23J), remarking that this approach has 
been also applied in the needlet framework to local thresholding by [3] and [J_2] . 

Section [2] will recall some preliminary notions, as needlets, their main prop- 
erties and the Besov spaces. Section [3J will describe the block thresholding 
procedure we suggest for needlet regression estimation , while Section |3] will 
present the main minimax results. Section [5] will collect some useful auxiliary 
results, while Section [5] will exploit the proof of the main result of this work, 
named as Theorem [T] 

2 Background results 

In this Section, we will review briefly a few of well-known features about the 
Voronoi cells on the sphere, the spherical needlet construction and the Besov 
spaces. 

For what concerns Voronoi cells, we are following strictly [3J: further de- 
tails can be found for instance in the textbook [34], see also [2] and [37]. From 
now on, given two positive sequences {aj} and {bj}, we write that aj ps bj if 
there exists a constant c > so that c~ 1 aj < bj < caj for all j. Furthermore, 
B XQ (a) — {x G § 2 : d(x,xo) < a} and B Xo (a) = {x € S 2 : d(x, xq) < a} de- 
note respectively standard open and closed balls on § 2 around xq S § 2 , while 
| A | is the spherical measure of a general subset A C § 2 . Given e > 0, the set 
S e = {xi, xn} of points on § 2 , such that for i ^ j we have d {xi, Xj) > e, is 
called a maximal e-net if it satisfies d (x, H e ) < e for x G § 2 , U Xi ^ F B Xi (e) = § 2 
and B Xi (e/2) (~l B Xj (e/2) = 0, for i ^ j. For all Xi £ S e , a family of Voronoi 
cells is defined as 

V (x^ = {i £ § 2 : for j 7^ i, d {x 1 Xi) < d (x, x^)} . (1) 

In [3J it is proved that: 

B Xi (|) cV(xi)cB Xi (e) . 

We resume now some features on the scalar needlet construction, system, 
suggesting for a more detailed discussion [3j5], [37J, see also [3] and [33]. A 
needlet system describes a well-localized tight frame on the sphere: it is a well- 
known fact (cfr. [36]) that any function belonging to L 2 (S 2 ) can be represented 
as a linear combination of the components of that frame, preserving furthermore 
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some fundamental properties of nccdlcts. Indeed, let us recall that the space 
Li (S 2 ) of square-integrable functions on the sphere can be decomposed as the 
direct sum of the spaces Hi of harmonic polynomials of degree I, spanned by 
spherical harmonics {Yi m y whose definition and properties can be found 

in [33] and [3] ■ If we consider 

l'-0 

the space of the restrictions to § 2 of the polynomials of degree less (and equal) 
to I, the following quadrature formula holds (see for instance [4]): given I G N, 
there exists a finite subset \i such that a positive real number (the cubature 
weight) corresponds to each £, £ Xi (the cubature point) and for all / G II;, 

Given B > 1 and a resolution level j , we call Xj B 2(j+i)j = Zj, card(Zj) — Nj] 

since now any element of the set of cubature points and weights, Xjk} , will 
be indexed by j, the resolution level, and k, the cardinality over j, belonging to 
Zj. Furthermore, we choose {ZjJ^^o be nested so that 

Nj &B 2 i,\ jk &B- 2 i . (2) 

We consider a symmetric, real-valued, not negative function &(•) (see again |4]) 
such that 

1. it has compact support on 

2. be C°° (R); 

3. the following unitary property holds for |^| > 1: 




For each ^ - fc £ Zj, given b (•) and B, scalar needlets are defined as: 

B3- 1 <(<BJ + 1 ^ ' 

where Li ( (x, y)) = Y, m =-l Ylm ( x "> Ylm (v) > describing loosely speaking the weighted 
convolution of the projection operator L\ {{x,y)). 

The properties of the function b (•) reflect three basic features of the needlets. 
Indeed, from the infinite differentiability of &(•), we have a quasi-exponential 
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localization property (see for instance [37]), which states that for k G N, there 
exists Ck such that for x £ § 2 

HM*)|<- 7^ 77T> (3) 

(l + B^^z))* 

where d (£j k ,x) is the geodesic distance on the sphere. In view of this property, 
it is possible to fix a bound (upper and lower), for the norms of needlets on 
L p (S 2 ), for 1 < p < +oo. Given p, there exist two positive constants c p and C p 
such that 

c p B^-l)<\\^ k \\ Lp{S2) <C p B^) . (4) 

Because the function b(-) has compact support in it follows that 

b (-gy) has compact support in \B$~ 1 , , hence needlets have compact sup- 

port in the harmonic domain. 

Finally, the unitary property leads to the following reconstruction formula 
(see again [36]): for / G L 2 (§ 2 ), in the L 2 sense, 

f(x) = Y,Pi^jkW , (5) 
j.fc 

Pjk '■= </i ^jfc) L2( g2) = ^ V'jfc («) / 0) cte , (6) 

where /3 j7s are the so-called needlet coefficients. 

Before concluding this Section, we recall the definition and some main prop- 
erties of the Besov spaces, referring again to [4], [12] and [23] for further theo- 
retical details and discussions. Let f £ (§ 2 ); we define 

G k (f,ir) = H MJf-H\\ Lw( p ) , 

which is the approximation error when replacing / by an element in . 
The Besov space S£ is therefore defined as the space of functions such that 
/ G L<x (§ 2 ) and 



oo \ 

£-(FG fe (/,7r)r <^ 

k=0 ) 



\k 

The last condition is equivalent to 

' oo \ 

Y / (B' r G B] (f,n)) q \ <oo. 
K i=o J 

Moreover, F £ B^ q if and only if, for every j = 1,2,... 

(£(|0ik|fell Mf p)) =£ i B 



-jr 
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where Sj G £ q and B > 1. The Besov norm is defined as follows: 



\\f\\ LAS2) +suM r+ ^ (P jk ) k 



q < oo 
q = oo 



As shown for instance in [4], if max(0, 1/n — 1/q) < r and n,q > 1, then we 
have 

f eBl q &\\f\\ B r 9 <oo . 

The Besov spaces present, among their properties, some embeddings which will 
be pivotal in our proofs below. As proven in [3] and [15], we have that, for 
tti < 7r 2 , q\ < q 2 



nr r nr W r K r E r C E " 1 1 " 2 



(7) 



3 Needlet Block Thresholding on the Sphere 

In this Section we will present needlet estimators for nonparametric regression 
problems and, then, we will suggest a procedure to fix blocks for any given 
resolution level j and, consequently, we will define the so-called needlet block 
threshold estimator. The first step is close to the one described in [3], [T2J 
for local thresholding, the other one being an adaptation to the sphere of the 
procedure developed for R in [24], [25], see also [23] 

In order to introduce the nonparametric regression estimator, let us initially 
define the so-called uncentered isonormal Gaussian process with mean / (see 
********). More precisely, we assume to have a family of Gaussian variables 
such that for all hi, /la € Sj, {X (hi), are jointly Gaussian with mean 

EX(h) = (h,f) = [ f{x)h(x)dx 
Js 2 

and covariance 

E (X(hi) - EX {hi)) (X(h 2 ) - EX(h 2 )) = (hi,h 2 ) L2(s2) . 

We shall in fact be concerned with sequences {X n } of such processes, where we 
assume that 



EX n (h) = (h,f) = / f{x)h(x)dx 



s-' 



and covariance 



, [X n (hi) - EX n (hi)) (X n (h 2 ) - EX n (h 2 )) - - (hi,h 2 ) L 2 {s2) 
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Consider now the usual needlet system {ipj k }j k an d let / € L P (S 2 ); we have 
the following: 

Pjk = EX n(ipj k ) = (>PjkJ) = f(x)xp jk (x)dx , 

Js 2 

Pjk = X n(lpjk) = Pjk+ £ ik;n , (8) 

where 

Ee jk;n = E {X n (iP jk ) - EX n (iP jk )) = , 

1 1 2 

Es jk;n = ~ (^jk^jk) L 2 {S 2 } = ~ \\lpjk lli2 (S 2) . (9) 

EEj kVin Ejk 2 ;n = — yPjki' ' t l } jk 2 /L 2 (S 2 ) 

In a formal sense, one could consider the Gaussian white noise measure on 
the sphere such that for all A, B C S 2 , we have 



EW(A)W(B) = [ dx 
J AnB 



so that 

£jk;n = ~ I ^jk( X ) W ( dx ) 



n 



s 2 



As described above (see also [4], [12]), / can be described in terms of needlet 
coefficients, up to a constant, as 

f = J2J2^^ 3 k ■ 

j>0 k=l 

Let us now define the blocks on which we will apply the thresholding pro- 
cedure: as anticipated in the Introduction, differently from [24], the structure 
itself of the needlet framework puts in evidence a way to build them. Let us fix 
j > 0: recall that for each resolution level j, we have Nj ~ B 2 ^ cubature points. 
Given the size of the blocks, i.e. the number of cubature points belonging to 
each of them - let us say £j - we will build using ([lj a set of Voronoi cells, 
containing lj cubature points. For each cell, we choose a cubature point £- a 
to index it: we define Sj (£j) as the number of Voronoi cells obtained to split 
cubature points into groups of cardinality £j . Let us define the set 



Rj < s = {k:Z jk eV(Z js )}, > l s 



j ■ 



From ([T]) , it is immediate to see that each cubature point belongs to a unique 
Voronoi cell. We choose £j such that 

t, = [NV] « N] . 
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where [•] denotes the integer part and < r\ < 1, such that 

N, 



5, = ^«(^y- 



3 

We finally define, for any integer p > 1, 



• ! -'- ; /* : ~p. E P P jk 



3 keRj- B 



and its corresponding estimator 



3 fee-Rj;, 



p 

jk J 



similar to the ones suggested in [51], Remark 4.7. 
Let us define the following weight function 



n I ' 



we have: 

/* = EE( E hkM «w > (io) 

where: 

• J n is the highest resolution level considered, taken such that 

B Jn = n' ; 

• k is the threshold constant (for more discussions see for instance [4], [12]): 

• the scaling factor t n , depends on the size of the sample. We will fix 

t n = vT^. 

4 Minimax ZArisk rates of convergence 

Our main purpose is to describe the performance of the procedure and the 
optimality of its convergence rates with respect to general L p (S' 2 )-loss functions. 
This result is achieved in the next theorem. 

Theorem 1 Let f 6 B r Kq (G), the Besov ball so that ||/||g r ^ < M < +oo, 

r — — > 0. Consider f * as defined by ilO\) . For p £ N, there exists a constant 
Cp — c p (p, r, q, M, B) such that 



Sup E ||/* -/Hip/ga^ < CpTl 



-a(r,7r,p) 
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whe 



rp 



. 21, - i : for ir > 2(j 2 ^ 2) 

a(r,n,p)={ p ( r - 2 (i-i)) 2p 
2( r _a(i-§)) /or ff < 5(?+2y 



If p = +00, i/jere existe a constant 

Coo — Cqo 

sup Elir-Zll^^coon-^'^ 
/es; (G) 



where 



a(r,7r,p) 



2(r-2(i-|)) 



Remark 2 T/iis procedure achieves minimax the rates established in J^j and 
fW, see also Ml- 



5 Auxiliary Results 

In order to prove Theorem [TJ we will need the following 

Lemma 3 Consider f3j k as described m[31 There exist constants CpjC^^C/ 
such that, for < , j = 0, J n , 



E 



(11) 



sup 

k=l,...,Nj 



Pjk ~ Pjk 



^CooU + ifn-p/^pyi , 



and for all 7 > there exists k > such that 



(12) 



A ■ — A ■ 



> «tS ) < Q 



(13) 



Proof. First of all, consider that, from Equations jSJ and ((9|), we have 



(l^-^f) = E(| £ife; „n 



(Var (ejfe ; „)) : 
1 



2fr(H±i) 



0(n~*) , 



I'' 



2fr(H±i; 
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to obtain (TTj!- Now, for Mill's inequality, if Z ~ N (0, 1), we have P (\Z\ > x) < 
^/2/7rcxp (— x 2 /2) jx , which leads us, from (|9j), to 

P(|ei*jn| >ar) = 2P(|Z|>vW) 

V 7r y/nx 
< C £ e-^ , 

On the other hand, we have 



E 



sup 

k=l,...,Nj 



x p_1 P 



sup 











Pi* 


> ay da; 


£jk;n | 


> x 


da; 



where 



£ 2 



c^dx , 



= C I x p ' x B^ : maxP(|£ jfc;n | > a;) da; 



31 ^ vTv 3 



We can easily see that 

while on the other hand, considering that for x > 2-%/2/nj 



we obtain 



E 2 C 



i?- 1 B v e- a *-dx 



so we achieve (|12j) . It remains to prove (|13j) . which corresponds to prove that 
for all 7 > there exists k > such that 



Pr 



1 3 
1 v^^p 



i/p 



< C„ 



1 



We start from rewriting 

/^■fc = VnPjk + Vnsjkyn = VnPjk + £jk 
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So we need to analyze terms of the form 

f £ * + ^ X>^r + - + ^ £ "T ■ <») 

* J fe=i ^ fe=i * J fc=i / 

We start by observing that: 

u \ 1/2 / . u \ 1/2 



J fc=l \ J fe=l / \ 3 k=l 



Now we know that 



fc=l k=l 

On the other hand, by Lemma 21 it is easy to see that for all p, 7 > 0, there 
exists k > such that 



It is immediate to see that 



by choosing suitable s and 7, we have 

— ^ — E ^fe = f E e .fc 

3 k=l \ 3 k=l 

The same holds for all the other mixed terms in Equation (fT4|) . ■ 
Lemma 4 Assume that EEjk = 0, Ee^ k = 1, and 

C 

Ee jki£jk 7 < -tt , for all M > . 

T/ien /or p € N, 7 > there exists k > smc/i £/ia£ 
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Proof. Without loss of generality we can take p to be even; note indeed that 



fc=i 



fe=i 



Now we can write 



fc=i 



T=l 



whence 



Pr I i E ^ > p( K + Ee %) h E Pr <! f E ^(^) > 



fe=i 



K 
Cr 



By a simple application of the Markov's inequality, the result will hence follow 
if we simply show that 



E 



k=\ 



< 



C 



Now let us take for notational simplicity r = 2; the argument for the other 
terms is identical. We have 



E 



1 3 



k=i 



1 

Jf E E { H -r( £ jki)—H T (Sjk y )} 

3 k\ j . . , k-y — 1 



1 \ 3 

= 77 S E I^fo*!^)]' 



7/2 



J-2 



+ 7r \ E [ S ( e J*i e jfei)] 

1 f ^ 

+7? ^ E [ £ ( e jfci £ ifei)] 



k\...ki 
-4 



J I fcifc 2 



X] E ( £ jki £ jk 2 )E(e J k 2 e J k 3 )E(e : jk 3 e j k i )E(e j k i e : jk 1 ) 

^ E { £ jk 1 £ jk 2 )— E { £ jk e £ jk 1 ) 
i. ..he 

if" 3 \ 

+ fl\ E E ( £ jki £ jk 2 )- E ( £ jk y £ jk 1 ) > 



= o(^7 7/2 ) + o(*~ 3_1 ) + ... + o(^ 7+1 ) , 
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because 



^ E { £ 3k l £jk 2 )--E{e j k^ej kl ) < E \ E ( £ jki£jk 2 )\ — \E(e jki _ 1 s jki )\ 



k\...k„ 



7-1 



6 Proof of Theorem [T] (upper bound) 

First of all, we stress that where explicitly mentioned, proofs will be totally 
similar to the standard thresholding case described in [3], hence we will omit 
them. Considering that JJs=i Sfce-R = 2k=i> ^ * s easv to see that 



'II/" /I1lp(s 2 ) 



j=0 s=l V k£R ji3 I j>0 fc=l 



Lp(S 2 



JV, 



= E 



< 2 P ~ 



j=0 g=l k£Rj- s 



E 



EE E (vAk-Pjk)i>jk 

j=0 s=l fcS-R J!s 



j>J„ fc=l 
P 



LP(S 2 ) 



Lp(S 2 ) 

EEm. 

j>,/„ fe=l 



L"(S 2 ), 



/ + // . 



CASE I: Regular Zone 

Consider p < +oo. For p < it, we have BJ g C £>p g : we can therefore take 

r-l + 2. 

tt = p . Consider instead the case p > it: we use the embedding B^ q C B pq p * , 
and moreover we assume 



2 r 



rp 



< 



p 2r + 2 2(r+l)p ~ 2p 



we have as in [3] that 



as claimed. 

Consider now the variance term. First of all, we easily obtain from the Loeve 
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inequality 



I < CE 



Jn Si 

EE E ( w ispjk- Pjkj^ik 

j=0 s =l k£R j;3 



Lp(§ 2 ) 



< cjf- 1 

3<J n 



E E 

s=l keRj, s 



Lp(S 2 ) 



As described in 0], see also [33], we have the following needlet property: 



E 



Hence we have 



E a ^. 



II^j*;|Ilp(s 2 ) 
Z>(§ 2 ) fc 



" E e k iip 



L?(S 2 ) 



J<Jn 



E E ( w 3>0jk- pjk)i> 



s=l fcS-R, 



Lp(S 2 



= E E 

3<Jn 



Si 



A 



+ E E 

3<J« 

+ E E 

J<Jn 



- E E 



E E ( w J'hk-Pik)w( 
E E ( w i>Pjk- Pn)i>iki\ 

s=l k£R j[s 

E E ( w 3'hk-^ik)^iki{ 

s=i feeflj-., 

E E ("v v v)'> ; ( 



.4 



Lp(S 2 ) 

p 



><£ I |A,.|< 



A; 



2 

<#JJ(|A ia |>2t£) 
<^)/(|A ia |<2*S) 



< 



c EE e imil^ (hk-Pik) i 



A 



3S 



Lp(S 2 ) 



Lp(S 2 ) 
P 



Lp{S 2 ) 



>tl)l(\A js \>^ 



j<J n s=l k£R j;s 
5 . 

EE E IWIL w E f(%*-^O Pj (W^*) J 0^'i < f 
EE E IIM^ 2 J/^r E k(^ 



y jfc|| i p(g2 ) E 

- E E E II^IIlp ( s 2 ) I^T e 

3<J„ s=l keRj [s 



A 



<^j/(|^ is |>2^) 
<t^J(|A js |<2^) 
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= Aa + Au + Ua + Uu . 

The idea now is to split this sum into four terms: in one of them, Aa, both 
the Aj s and Aj S are large; in another one, Uu, they are both small and in the 
last two of them, Au and U a, the distance between Aj S and Aj S is large. In 
the first two cases, in order to achieve the minimax rate of convergence, we will 
split the sum into two parts and we will show the convergence of each term by 
using mainly (j4|), (fTT|) and (|12[) . The convergence of the last two terms will be 
proved by using (|13|) . We start by studying 



j<J n s=l k£Rj. s 

< ^EE B^-Vl(\A js \>% )E 

j<J n s=l keRj; S 



Pjk — Pjk 



As in [3], we fix J\ n such that 



hence we have 



Aa < Cn-vl* £ ^ , /r , */ (|^,| > |) + £ ]>>i^- 2 >J (|A, S | > | 



v., 



\j<Jm j=Jm «=1 





' 11 


l> 






2 




)) 
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Observe now that 

E jrt jB ^-vi (\A js \ > |) < E tjB"*-* £ 1^.1 (| 

j = J ln 8=1 ^ ' j=J In 8=1 V 



< ct-j £ e jB ^J2l e 

j=Jln 8=1 J fcei? 3 - 

< a-p e s^^Ei^r 

3=Jl» fc=l 

< cn« e Ei^r^" 2) - 

j=Ji„ fei=i 

Now, because / 6 £?p 9 , we have 

at, JV 3 

e i^r^- 2 * = c E i^r im* - cB ' pr3 ■ 
fc=i ^=1 

Hence, as in [4] 

n5 E E Bj(P ~ 2) < Cn^TTJ < gpJm . 

j=Ji„ fci=i 

Finally 

Consider now the term Uu. We have 

Si 

Uu < CEE E fell^) 1^*1^(1^.1 < 2^) 

j'<J„ 6 = 1 k£R j:s 

< cE^ (p - 2) E E l^r^(i^i<20 

j<J„ S = lfce-Rj; S 



s=l 



< c 



C 



E NjB^-^ti + e ^' (P ~ 2) E^ 



E^' (p - 2) -- f +EEI^*r ll^ll^) 
<j ln j=j ln k=i 



n 2 ('-+ 1 ) 
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Let us study now Au and Ua. As in [3], we have 



LJ>(§ 2 ) 



^ * E E E fell* 

j<J„ s=l k£R j]s 

< CB pJ "n-%n-' 1 < Cn-~< 
Si 



E 



2/; 



4 ■ — 4 ■ 



> 



ua < E E E Ik 

Because for r > 1, we have 



J^IIlp(S 2 ) 



71 7 < 71 2 < n 2(r + l) 



the result is proved. 

Consider now p = +oo: we assume now / £ ^j, to obtain 



E|ir-/IL<: 

As in [1], we have: 



j=Q k=l 



L°°(S 2 ) 



A', 



E E^^'fc 
3>j n fe=i 



= /+!/ . 



L°°(S 2 ) 



II = o(n 2 <'+ 1 > 
For what concerns /, we instead have 



A', 



E ( W ^3k ~ ^fe) ^fc 



fc=l 



< 



L°°(S 2 ) 



J=0 



sup [WjP jk - /3 ifc 



< c£b'E sup(^ fe -/3 

J» r , 

su p - Pik) A 

3=0 1 k V 

Jn 

+C E sup |/3 jfc |E[/( - A, 



KJl 2 



> 



KI1 2 



3=0 
Jn 



> ku 2 



+CE SJ sup| i 3 jfc |/(|A 7 -| < 2ktT 
j=0 fc 
= Aa + Au + Ua + Uu . 

Again, we choose Ji jn such that 

B Jl - n = k'ti^ 1 ) ; I ( \Aj\ > 



for j > Ji „ , 
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and similarly to [2], we obtain 

Aa < CJ hn n~^B Jl " < Cn~^+v ; 

The other two terms Au and U a are similar to the case previously described. 
For general it and q, we observe that B r ^ q C B^^, r' — r — 2/n. Hence we 
obtain 

W - /lli=o C s2) < C^rT^ 1 ) = CJ n n-^-wt-im> , 

as claimed. 

CASE II: Sparse Zone 

The proof follows strictly the same steps of the regular case. Indeed we have 
B^ c S P9 . Hence we have 



■||/*-/llL>(s») < 2"-'E 



j=o s=i fce_R ; 



Lp(S 2 ) 



A', 



E E^ 



= :I + II. 

Also in this case, as in [4] , because r — 2/ir + I > 1, we have for the bias term: 
II = 0(n-^-k)V^-^-k))) . 

On the other hand, we split / again into four terms as above. On one hand, we 
obtain 



Au ^ EE E IMIW) E 



j<j„ s=i keRj 

Si 



fijk ~ Pjk 



J 2p 



Aj Aj 



> 



ua < EE E II^IIW)l^r( p ( 



j<J„ s=l keRj 



whose upper bounds recall exactly the same procedure developed in regular 
zone. On the other hand, consider initially: 



Aa < CnS Y, E E B^l(\A js \> 

j<J n s=l keRj ;s 



In this case, we fix Ji n so that 



for j > J 2n ■ 
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to obtain 

Aa < Cn-i E E E £ j(p - 2) /(Wl>f) (15) 



J<J2n 



Ni 



< en-**- 1 x: £ j(p - 2) <MEl^l 



j<J2n \fe = l 

' Ni 



j<^2„ \fe = l 



AT,- 



3<Jln 

*+'(*-*)) 

because in the last inequality we used 

l- p (p-2-(r+l-*))_ f(r + 2(i-i 



2 2(r-| + l) 2(r-§ + l) 

Consider now 

£/u < C ]T B*- 2 )^4 is /(|A j8 | < 2i£) = C ^ B^'^Vd^l < 2**) 

J<J„ S=l 3<Jln s=1 

+ g sid'-a)^ £ A js / (|A,- a | < 2tl) = Uux + Uu 2 . 

j = J2n 8=1 

Now fix 

P-2 

m = 2 7 i 

r - | + 1 

so that 

p-m = p V— — >0; 

r - t + 1 

p-7r(r + l) 
m-7r = 5 - >0 . 
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Simple calculations lead to 



Um = C Y, B^- 2 h j J2A js I(\A js \<2t* n ) 

j<J 2n 8=1 V tn 

Sj j 

J<Ja» 8=1 

< ct^ Yl B^tpYl E r» 

j<j 2n s=i yfceRj a 

j<J2n \fe=l 

< a- 1 E ^ 2) ^fEi^rii^iir^)) 

< CtP- 1 B j ( p - 1+ ™) B- jr 

j<J2n 

= CtP- 1 J2 B^p- 2 -^- 1 -^)) 

j<J2n 

= ctp-^M?- 2 -^- 1 -!)) 

< Cn^B^- 2 -^- 1 -^)) , 

exactly as above in Equation (fT5j) . We have to study just the last term 

Uu 2 = C Y NjB^P-^AjIHAjl < 2t p n ) . 
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Wc have 



Uu 2 = C £ B^l j Y,A js I{\A js \<2tl) 

j = J2n 8=1 
Jn 

< c E ^ (p " 2) ^E^/.V PJ (i^i< 2 ^) 

j=J2n 8=1 

< c*jr™ E Bf- 2) ^X>g 

J = J2n 8=1 

< ctr™ £ ^ a) 4 _ 'i;( E l^r 

j = J 2 „ «=1 \k€R js 

Jn ( N 3 

< ai- m lis, Ei^r 

j=j 2n \fe=i 

< c*r m E BJ(P ~ 2) & • ^) (e l^rl 

j=J 2n V fe / 

< en- £ ^- m) fEi^ni^nr m( sO 

j=J 2n V fe / 

< C't P ~ m ^ #7(P _m )_B _ -? m ( r_ f +m) 

J = J2„ 

We can easily see that 

( 2 2 \ p- 2 / 2\ 
(p-m)-mr-- + — = p - 2 ^ 1 + r - - =0 

Hence 

K^(l-i)) - 
f7u 2 < Ct£-" = O ( n 2 ("-* +1 ) 



as claimed. 
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